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Abstract 



We present a new software tool to enable astronomers to easily compare observations of emission line ra- 
tios with those determined by photoionization and shock models, ITERA, the IDL Tool for Emission-line 
Ratio Analysis. This tool can plot ratios of emission lines predicted by models and allows for compar- 
ison of observed line ratios against grids of these models selected from model libraries associated with 
the tool. We provide details of the libraries of standard photoionization and shock models available 
with ITERA, and, in addition, present three example emission line ratio diagrams covering a range of 
wavelengths to demonstrate the capabilities of ITERA. ITERA, and associated libraries, is available from 
http : //www . brent groves . net/itera . html. 
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1. Introduction 

All emission lines are sensitive to some extent on 
the physical conditions in the emitting gas. This in- 
cludes the density and electron temperature of the 
medium, the abundance of the emitting element, 
and the details of the ionizing source such as the 
luminosity and effective temperature of the star in 
Hn regions. While both atomic and line theory 
and observational data have advanced sufficiently 
over the years such that the dependency of individ- 
ual emission lines on nebulae parameters is known, 
determining these same parameters from individ- 
ual lines is difficult apart from a few exceptions. 
One of the simplest ways to break these inherent 
degeneracies is through the use of emission line ra- 
tios. Using lines arising from the same elements, or 
similar ionization potentials, or even from different 
levels of the same ion, can minimize the dependence 
on some or most of the parameters that control line 
emission, thus enabling the construction of diagnos- 
tics. 

One of the first diagnostic emission line ratios 



was the 
ward by 



O III] line rat i o A43 63A/A5007A put for- 
Menzel et al. (|l94ll) . By using forbidden 



lines of the same ion, +2 , both the abundance 
and ionization dependencies are removed, leaving a 
ratio strongly dependent on electron temperature 
due to the different excitatio n levels of the lines in- 
volve d (see e.g. Fig. 2.5 in iDopita fc Sutherland! 
2003). This line method, and other methods of 



determining electr on temperature , have been dis- 
cussed in detail bv lPeimbertl(|l967t ). The same rea- 
soning has also been extended to several other line 
ratios, using lines from a single ion but arising from 
different energy levels, such as optical-IR ratios 



like [ O in] A5007A/A52/xm ( see e.g.lDinerstein et al 
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119851 ) and UV ratios like C in] A1909/C iiiA977. 

Similarly, ratios of forbidden lines arising from 
the same ion and similar energy levels remove abun- 
dance, ionization and temperature dependencies. 
If these lines arise from levels with different col- 
lisional de-excitation rates or radiative transition 
probabilities, this ratio will be sensitive to the den- 
sities in the range spanned by the critical densi- 
ties of the l ines involved. Th is idea was first put 
forward by lAller et all (|l949l ) for the [On] dou- 
blet, A3726/A3729, and has been extended to sev- 
eral other strong lines like [S n] A671 6/A673 1 or in 
the IR [S in] A18.7/im/A33.5/mi (see iRubinl . Il989l 
for a nice overview of these ratios) . 

Apart from the temperature and density sen- 
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sitive ratios described above, determining physi- 
cal quantities from single line ratios proves dif- 
ficult due to degeneracies between the various 
controlling parameters of the nebula emission. 
Generally what is done is to use several line 
ratios in combination to break these degenera- 
cies, such as determining the m etallicity or abun- 
danc es in the gas (see e.g. Kewlev fc EllisonL 
120081 . for an overview of several line methods for 
metallicity determination), or to use empirically- 
based assumptions, such as the constancy of the 
Balmer decrement in H 11 regions for dust-reddening 
dete rmination (also t heoret i cally supporte d, see 



Cox fc Mathewsl. Il969i: iBrocklehurstl Il971 



iMiller fc MathewsT 19721 though these concentrate 



mostly on planetary nebulae). An overview of all 
methods of comparing ob served emi s sion li nes with 
theory can be found in IStasihska ( 2007) a nd in 



books such as Osterbrock fc Ferlandl ( 2006) and 



iDopita fe Sutherland! ( 2003 ). which give a broad 



understanding of the topic and underlying physics. 

One of the best ways to disentangle these de- 
generacies is to plot two line ratios against each 
other in what is commonly called a line diagnos- 
tic diagram. These were first used as a method to 
di stinguish various clas ses of emission- l ine ga laxies 



Heckmanl (|l980t) and lBaldwin et all (|l98lh . with 



the diagrams described in the latter still commonly 
used to distinguish star-forming or H II galaxies 
from galaxies dominated by an active galactic 
nucleus (AGN). These diagrams for emission- line 
galaxy classification w e re rev ised and extended by 



Veilleux fc Osterbrockl ([19871 ). who based their se- 



lection of line ratios on 5 criteria: Line strength (1), 
Line proximity (for blending (2) and reddening (3) 
effects), the preferable inclusion of hydrogen lines 
(4) and the capability of detecting these lines (5). 

With improvements in the availability of spectra, 
with wavelengths from the UV (such as Hubble- 
STIS) to FIR (e.g. Herschel -PACS & -HIFI) avail- 
able, increased spectral sensitivity and resolution, 
and a copious number of ion ization models now 
available (e.g. iMorisseti 120081 ). such criteria need 
not be so strict. However, these large increases 
in wavelength range, sensitivity, and models make 
analysis of emission-line objects much more compli- 
cated. 

It is for these reasons that we introduce here a 
new, publicly available, emission-line analysis tool, 
ITERA. ITERA, the IDL Tool for Emission-line 
Ratio Analysis, allows the simple comparison of 
observed emission-lines from UV to far-IR wave- 



lengths to models to assist in the determination of 
properties such as gas density, metallicity and ex- 
citation mechanism. We describe in the following 
sections the code and demonstrate some of the pos- 
sible diagnostic diagrams available. 



2. The ITERA Program 

The IDL Tool for Emission-line Ratio Analysis 
(ITERA) had its inception with the SH OCKPLOT 
tool cr eated to display the results of the Allen et"al1 
(|2008f ) 150-1000 kms^ 1 shock models. The power 
of such a tool was quickly realised and the original 
tool was rewritten and expanded into its current 
form. ITERA is an IDI0 widget tool which enables 
astronomers to plot line ratio diagrams of emission 
lines arising from atomic and ionized species as de- 
termined by standard photoionization and shock 
models. It also allows for a comparison of these 
line ratios with those obtained from observations 
or other models. As discussed in the introduction, 
these diagrams, or the individual ratios of the di- 
agram, can be used as diagnostics of nebulae exci- 
tation mechanisms, nebulae density, electron tem- 
perature, metallicity or individual abundance vari- 
ations, and other model-specific parameters. IT- 
ERA can also be used to determine line sensitivi- 
ties to such parameters, thus giving the possibility 
of creating new or "best possible" diagnostics for a 
given set of observations. It can even be used to 
estimate line fluxes that are weak or in unobserved 
wavelength ranges. 

ITERA works in three simple steps; the user se- 
lects the emission lines to be used as the denomi- 
nators and numerators of the ratios for the X- and 
Y- axes, selects the model(s) that the user wants 
to examine, and, if wished, the user enters obser- 
vational data for comparison. These steps can be 
performed in any order. Multiple lines can be cho- 
sen for the ratios in the case of emission line dou- 
blets or line blends due to resolution or velocity ef- 
fects. What is returned is an emission-line ratio di- 
agram of the chosen ratios, displaying a grid of the 
chosen models and, if entered, observational data 
points. This diagram can be printed, or the model 



^DL, the Interactive Data Language, is a 
computing environment for data analysis, vi- 
sualization, and application development, avail- 
able from ITT Visual Information Solutions 
(http: / /www. ittvis . com/ProductServices/IDL. aspx ). 
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grids displayed can be output for use in other rou- 
tines, or direct comparison with data. The program 
and associated libraries and routines are available 



from http : / / www . brentgroves . net/itera . html 

An IDL virtual machinqj version is also available 
for those without access to IDL licences. Full in- 
structions and details of the program are found on 
the webpage, but we describe below the major parts 
of the program. 

The principle of ITERA is to allow ease of ac- 
cess for astronomers to existing shock- and photo- 
ionization models for the determination of emission 
line strengths. Thus the strength of ITERA lies in 
the library of shock- and photo- ionization models 
associated with the code. This library is a distinct 
yet necessary part of ITERA, needed to generate 
the line ratios. As the library is not an integral 
part of ITERA, only the desired model series need 
to be downloaded and installed (such as AGN or 
Hn-region models, with further models added later 
as required, including photoionization models gen- 
erated by the user. 

2.1. Choosing the Models 

The library of models associated with ITERA 
cover the range of excitation mechanisms expected 
within emission-line galaxies; star formation, active 
galactic nuclei and shocks. These models, gener- 
ated by the shock and photoionization code MAP- 
PINGS in, broadly cover the ratio space occupied 
by the dominant optical emission lines from galax- 
ies. Other sources of emission lines, such as individ- 
ual O and B stars (Hn regions) or planetary neb- 
ulae, can be included through the user-generated 
option of ITERA, and may be included in the li- 
brary at a later date. The ITERA model library is 
split into four main categories: Starburst or stellar 
ionization models, AGN models, shock models, and 
user-generated models. The currently available li- 
braries are shown in table [TJ and described briefly 
below. 

ITERA has the possibility of plotting up to any 
number of shock or photoionization model grids, al- 
though for clarity the number displayed on any plot 
is generally limited to 5. Every model set has a lim- 
ited number of parameters, with common parame- 
ters to the models being the emitting gas metal- 
licity or abundance pattern, and the gas density 



or pressure. The other parameters associated with 
the models tend to be specific to the specific sets, 
either due to the different physics associated with 
the ionizing/excitation mechanism, or due to differ- 
ent definitions by the model set creators. These in- 
clude the hardness of the radiation field, the veloc- 
ity of a shock ionizing the gas, or the dimensionless 
photoionization parameter. The parameters asso- 
ciated with the different model sets within the IT- 
ERA model library are described in the following 
subsections. 

The model grids displayed within the ITERA 
plot window are formed by the ranges of two of 
these parameters, with the other parameters of the 
model sets fixed. The fixed parameters thus de- 
scribe the different model grids that can be plot- 
ted. The first parameter of the model grid is set 
and cannot be chosen by the user, and is generally 
a measure of the ionization state of the gas, such 
as ionization parameter or shock velocity. The sec- 
ond parameter of the grid can be chosen by the 
user, with the choice of density, metallicity or the 
unique parameter associated with that model set 
(as described below). If the user of ITERA wants 
to compare parameters other than the set first pa- 
rameter, e.g the ionization parameter, this can be 
done by simply limiting the parameter range to a 
single value and plotting multiple grids. 

2.1.1. Starburst Models 

The ITERA model library includes 3 different 
Starburst models, aimed specifically at galaxies 
whose line emission is dominated by gas ionized by 
recent (< 10 7 yr) star formation, generally classified 
as starburst, Hll, or star-forming galaxies. How- 
ever, these models can also be generally applied to 
any massive star formation region, though the com- 
parison may break down for less massive star forma- 
tion regions once the effects of the stochastic sam- 
pling of the stellar initial mass function becomes 
important (< 10 4 M Q ). 



The first model set is from iKewlev et al.l (|200lh . 
These models use stellar ionizing spectra gen- 
erated by t he stellar population synthesis c odes 
PEGASE-2 dFioc fc Rocca-Volmerangel Il997t ) and 



2 IDL virtual machine 
time version of IDL 



Starburst99 (jheitherer et all Il999h . and consid- 
ers both continuous star formation and instan- 
taneous bursts. These models are then passed 
through the MAPPINGS in code for a range 
of ionization parameters^). The stellar models 



is free run- 
available from 

http://www.ittvis.com/ProductServices/IDL/VirtualMachine.aspjfThe ionization parameter, q, also known as the ioniza- 
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Starburst models 



Kewlev et al. (2001) models 

• Pegase2 stellar models 

o Instantaneous Star Formation 
o Continuous Star Formation 

• Starbust99 stellar models 

o Instantaneous Star Formation 
o Continuous Star Formation 
Levesaue et al. (2010) models 

• Standard mass loss stellar models 
o Instantaneous Star Formation 

o Continuous Star Formation 

• High mass loss stellar models 

o Instantaneous Star Formation 
o Continuous Star Formation 
Dopita et al. (2006) 7Z parameter models 



AGN models 



Standard, constant density, dust-free, 

photoionizcd AGN 
Dusty, Radiation-Pressure dominated, 

photoionizcd AGN 



Shock Models 



Shock ionized cloud 
Photoionizcd Precursor 
Shock+Precursor 



User-generated models 



User run Cloudy & MAPPINGS in 
photoionization models 



Table 1: Listing of the model libraries currently available 
with ITERA 



and subsequent photoionization modelling are de- 
termined for 5 metallicities, ranging from 0.05 
to 2.0 Solar (Zq). For full details see the 



Kewlev et al.l (|200lf ) paper. These models were 
used as a b asis for the strong - line m etallicity es- 
timators in iKewlev fe Dopital (|2002t ) and for the 



empirical / theoretic al AGN/star-forma t ion div iding 



lines discussed in iKauffmann et al. ([2003) and 
IKewlev et all (|200 " 



The second set of models is from iLevesaue et al 
(120101). and is an updated version of the 



Kewlev et al.l (|200ll ) models, using a later version of 



MAPPINGS in (version mr) code for the radiative 
transfer, and using the latest vers i on of the Star- 
burst99 code ( Vazquez fc Leitherer t 20051). L i ke the 



Kewlev et all (|200lf ) work. ILevesaue et all (|2010f ) 



model includes both instantaneous burst and con- 
tinuous star formation stellar models, and covers a 
similar range in ionizat ion parameter, q . and met al- 
licity. In addition, the Levesaue et al. ( 2010h mod- 
els also cover a range of ages between to 10 Myr 
(in steps of 0.5 Myr). What distinguishes these 
models from previous works is that they examined 
two different evolutionary tracks for the stellar pop- 
ulations, one with a standard mass loss treatment 
(STD) and one with high mass loss rates (HIGH) 
meant to more accur ately match observations (see 
ILevesaue et al. , 2010l for full details about the dif- 
ferent tracks and the models). 

The third s et is the "1Z parameter" models of 
Dopita et al. ( 20061) . As with the previous set, 
these use MAPPINGS nir for the radiative trans- 
fer and the latest version of Starburst99 to gener- 
ate the ionizing stellar spectra, though only instan- 
taneous bursts were considered. The same range 



of ste llar ages and metallicities as ILevesaue et al 



(l2010h are used. What differentiates these models 
from the previous sets is that a direct link with the 
stellar populations and the surrounding ionized gas 
is considered. In these models all stellar clusters 
are born buried within their molecular birth clouds, 
but, over time, they clear an empty bubble around 
them, sweeping up the gas that forms the surround- 
ing H II regions. This means that there is a direct 
link of both the gas density and incident ionizing 
flux, and thus the ionization parameter, at the in- 
ner edge of the Hn region to the evolution of the 



tion front speed, is a measure of the ionizing photon flux 
over the (hydrogen) gas density; q = S*/riH, where S* is the 
integral of the photon flux above 13.6 eV. This is also often 
given in a dimensionless form; U = q/c. 
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stars. With this link, the time-averaged ionization 
parameter is found to be directly correlated with 
the mass of the ionizing cluster divide d by the aver 



age IS M pressure, defined within the iDopita et al 



(2006) work as K (= M c i/P ISM )- Thus it is this 
1Z parameter rather than the local ionization pa- 
rameter, q, that is measured when the emission-line 
ratios are measured on galactic scales. For full de- 
tails about the 1Z parameter and the models see 
Dopita et all (|2006h . 



One important caveat with these models is that, 
while all three use the same global metallicities, the 
elemental abundances at each metallicity for each 
model differ due to the evolving definition of the 
abundance pattern at solar metallicity. Thus such 
differences must be considered when comparing 
these models. In table[5]we list the solar abundance 
patterns used in all models of the ITERA library 
discussed here. Most abundances are the same, 
except for oxygen which wa s recalibrated betwee n 
different model sets (see e.g. lAsplund et al" . 20051) . 
Different metallicities simply scale as multiples of 
the abundances listed here, except f or helium and 



nitrog en, a secondary element (see iGroves et al 



2004aj). One cautionary note on the abundances 



listed here is that these do not take account of 
the effects of dust depletion, which varies amongst 
the models (from non-existent to significant), which 
will alter the gas-phase abundances. 
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Element 


Kcwlcy2000 


Lcvcsquc2010 


1Z param 


AGN models 


Shocks 


H 


0.00 


1.00 


0.00 


0.00 


0.00 


He 


1.01 


-1.01 


-1.01 


-0.99 


-1.01 


C 


-3.44 


-3.59 


-3.59 


-3.61 


-3.44 


N 


-3.95 


-4.22 


-4.22 


-4.20 


-3.95 





-3.07 


-3.34 


-3.34 


-3.31 


-3.07 


Ne 


-3.91 


-3.91 


-3.91 


-3.92 


-3.91 


Na 




-5.75 




-5.68 




Mg 


-4.42 


-4.47 


-4.47 


-4.42 


-4.42 


Al 




-5.61 




-5.51 


-5.53 


Si 


-4.45 


-4.49 


-4.49 


-4.49 


-4.45 


S 


-4.79 


-4.79 


-4.79 


-4.80 


-4.79 


CI 




-6.40 




-6.72 




Ar 


-5.44 


-5.20 


-5.20 


-5.60 


-5.44 


Ca 


-5.64 


-5.64 


-5.64 


-5.65 


-5.88 


Fe 


-4.33 


-4.55 


-4.55 


-4.54 


-4.63 


Ni 




-5.68 




-5.75 





Table 2: Solar abundances used in ITERA Libraries, given as logio(X/H). 
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2.1.2. AGN Models 

The library of AGN models used in the ITERA 
code are those from the Groves et al. I (l2004allbh oho- 
toionization models of the narrow line regions of 
AGN. We include two series of AGN models, the 
classic dust-free, AGN power-law ionized models 
and the dusty, radiation pressure dominated mod- 
els introduced by Groves et al.. Both of these mod- 
els use simple power-law radiation fields (f v oc 
to represent the AGN ionizing spectrum, and ex- 
plore the same range in metallicities (with the same 
abundance patterns), gas densities and ionization 
parameters. However the first models are dust free 
and assume a constant gas density throughout the 
narrow line region cloud. The second set are taken 
to be isobaric (constant pressure) and include dust 
and the effects of radiation pressure on dust, which 
causes significant differences in the resulting emis- 
sion line spectrum at hig h ionization parameters 
(see iGroves et al. . 2004alfbl . for full details of the 
models) . These latter models are much more closely 
matched with the observed line ratio space and are 
the recommended model to use. 

While these models do not represent the full 
range of theoretical models for the emission 
from the narrow l i ne regions of AGN (see e.g., 
Binette et all . Il996t iFereuson et all Il997h , they do 
broadly cover the expected emission of strong emis- 
sion lines, and can be used to estimate the "AGN 
space" of line ratios. 



2.1.3. Shock Models 

The shock models used i n the ITERA libra ry are 
directly imported from the Allen et "all (2008) work 
on fast shocks. They include a range of metal- 
licities, densities, and magnetic parameters (the 
"unique parameter" of these models), and explore 
velocities from ~ 150 to 1000 kms -1 , where strong 
J-shocks occur. The shock models arc split into 
three sets. The first is the emission from the post- 
shock region, where gas is shock-excited to high 
temperatures and ionization states. The second 
model set is from the pre-shock or precursor region 
which is photoionized by the radiation emitted up- 
stream from the post-shock region. The third set 
is the combination of these (shock -(-precursor) for 
integrated measurements where the individual re- 
gions cannot be distinguished. For full details on 
the pa rameters and unde rlying physics of the mod- 
els see I Allen etaD (|2008l ). 



2.1.4- User Generated Models 

While the above models cover the expected 
parameter space of emission-line galaxies, they 
are not representative of individual objects, where 
variations in the abundances of specific elements, 
ionizing source variations or even geometrical 
effects can all affect the final emission-line spec- 
trum. Therefore we include the final option 
of user generated model libraries. This allows 
the possibility of a user inputting the results 
from well tested and resp ected photoio ni zation 



codes such as Cloud\0 dFerland et al 
and MAPPINGS n{f (|Groves et all 



19981) 



2004aD 



or using on-line interfaces such as mapiiionline: 
http : //www. if a. hawaii . edu/~kewley/Mappings/. 
For a listing, comparison, and bench marking of 
such models see lPeauignot et al. I (l200lh . though all 
these models have advanced since this work. 

To assist in the input of user generated models 
into ITERA, we provide tools for converting the 
outputs from these models to format used by IT- 
ERA. Currently only conversion tools for the pub- 
licly available codes Cloudy and MAPPINGS in are 
provided, but this may be expanded at a later date 
depending on demand. 

2.2. Choosing Line Ratios 

In principle, the choice of emission lines to use 
in the diagnostic line ratios in ITERA is limited 
purely by those provided by the code used to gen- 
erate the photoionization/shock model. In MAP- 
PINGS in (including all the models provided in the 
library) this is limited to ^1770 emission lines in the 
infrared, optical, and ultraviolet, and even includ- 
ing high-ionization/inner-shell lines in the X-ray. In 
Cloudy, this is further expanded by a host of weaker 
lines as well as molecular ro-vibrational emission 
lines from CO and H2. The lines are listed in sim- 
ple selection boxes, with only strong lines listed by 
default, though all lines, or lines in certain wave- 
length ranges, can be displayed if wished. 

In addition, sums of lines can be used in the line 
ratios to account for line blending at low spectral 
resolution, or to account for different levels or ion- 
ization states (i.e. [On]A3726 + 3729). 



4 Cloudy is available for download from 
http://www.nublado.org, along with full details on 
the code, how to use it, and how it works. 

5 MAPPINGS III is available to download from 
http://www.brentgroves.net/rnapiii.htinl, which in- 
cludes a description of the code and its use. 
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In practice however, the guidelines discussed in 
the introduction form the bases for good diagnos- 
tics. Lines of the same element and of the same 
ionization state or of similar excitation potential 
in different elements provide some of the strongest 
diagnostics. Lines relative to hydrogen recombina- 
tion lines are also good diagn ostics. As discussed in 
Veilleux fc Osterbrock (1987). lines in close wave- 



length proximity limit the possible instrumental, 
calibration, and reddening effects. 

Hence it is for this reason that several stan- 
dard, strong-line, diagnostic ratios, such as 
[S n]AA6716, 31/Ha versus [O m]A5007/H/3, are also 
included in ITERA as predefined diagnostic dia- 
grams for ease of use. 

2.3. Comparison with Observations 

Ultimately, the use of ionization models lies in 
their diagnostic power for observations, which is 
why ITERA also includes the ability to enter data. 
To allow for easy comparison of observations with 
the models, the ITERA program has three sim- 
ple ways in which data can be entered into the 
emission-line diagnostic diagrams directly; direct 
point by point entry, a simple list of line ratios for 
multiple data, and via tabulated emission line data 
for multiple line ratio comparisons. For the first two 
options, the line ratios for the chosen diagram are 
entered in and compared directly with the model 
on the chosen line ratio diagram (as shown in figure 
[2]) . Once the chosen diagram has cleared or another 
line ratio diagram chosen, the data are cleared. The 
tabulated emission line data are chosen in a similar 
manner as the models, with the data plotted over 
the model grids if the contain the chosen emission 
lines. Templates for the tabulated data are avail- 
able with the ITERA program. 

Conversely, the fluxes (relative to H/3) from the 
selected model grids for the four emission lines of 
the chosen line diagnostic diagram can be exported 
to formatted ASCII files for comparison with ob- 
servational data. This is especially useful in cases 
of large datasets, such as with SDSS. 

3. Emission Line Ratio Diagnostics 

To demonstrate the capabilities of ITERA, we 
present here three emission-line ratio diagrams cov- 
ering a range of wavelengths and diagnostic capa- 
bilities. These diagrams are all direct output from 
the ITERA routine, excluding the histogram in the 
first diagram. 




-2.0 -1.5 -1.0 -0.5 

log,.( Nil 6548.0A.6583.3A/ HI 6562.! 



Figure 1: The BPT diagram of [Nil] /Ha versus [Om]/H/3, 
showing grids of photoionization mod els for star form- 
ing galaxies (left. [ Levesque et all I2010T ) and AGN (right, 
IGroves et all teOCMbl') , described fully in the text. In both 
grids, lines of constant metallicity are shown by the blue 
curves, increasing from dark blue to light blue, and lines 
of constant ionization parameter are shown by the orange 
curves, increasing in value from red to yellow. In the back- 
ground is an image of a 2d density diagram from the SDSS 
DR4 sample described in lGroves et alj d2006f) . showing star- 
forming galaxies (green) , Seyfert galaxies (blue) and LINER 
galaxies (red). 



The first diagram (figure [lj shows [Nil] /Ha ver- 
sus [Qui] /H/3 the c lassic "BPT" diagram from 
Baldwin et al. (Il98lh that is a diagnostic of the 
excitation mechanism in emission-line galaxies. 
Shown on this plot are two model grids; the 



Levesaue et al. I (1201(1) STD models for star - forming 



galaxies on the left, and the IGroves et al.l (|2004bl ) 
models for Dusty AGN on the right. Both grids 
show a range of metallicities (blue curves) and ion- 
ization parameters (orange curves), though note 
that these grids are not equivalent. The exact pa- 
rameters that the star forming models have are: 
Z ~ 0.05, 0.2, 0.4, 1.0, 2.OZ and q = 1 x 10 7 ,2 x 
10 7 ,4 x 10 7 ,8 x 10 7 ,lx 10 8 ,2 x 10 8 ,4 x lO^ms -1 , 
for a 2 Myr old burst with density nn = 100 cm -3 , 
while the AGN models have: Z ~ 1.0, 2.0, 4.OZ 
with logU = 0.0, -0.3, -0.6,-1.0, -3.6, -4.0, a 
gas density of nn = 1000 cm -3 and are ionized by a 
power-law of index /3 = —1.4. Underlying the grids 
is a coloured 2D density histogram of the SDSS 
DR4 e mission line galaxy sample from lGroves et &T 



(2006), displaying the observed distribution of star 
forming galaxies (Green squares), Seyfert galaxies 
(blue squares), and LINER galaxies (red squares). 

The model grids clearly encompass the observed 
line ratio space. Note how the star forming mod- 
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Figure 2: The mid IR emission line ratio di- 
agram of [Ne m]15.6/^m/[Ne n] 12.8/^m versus 
[Slv110.5/jm/[N e n112,8 /jm, showing the same 

ILevesgue et al, I d201Cfl star forming galaxy photoion- 
ization model as in figurt^ with the grids curves as labelled 
in the colour bar to the right. The stars are a sample 
of O and B star ionized systems (star-forming galaxies 
and extraga lactic and ga l actic Hll systems) collated and 
described in Groves ct al. (2008). 



els follow closely the green branch of this diagram, 
showing the known strong dependence on metal- 
licity of the [Nil] /Ha ratio, with the ionization pa- 
rameter affecting the position of the models in both 
ratios, explaining the spread in observed galaxies 
(along with the average age of the ionizing stars). 
Also note that the models fail to reproduce the ob- 
served line ratios at very low metallicities, and do 
so even when different ages are used for the ionizing 
stars, indicating that there are still underlying is- 
sues with either the photoionization models, stellar 
population models, or both. The AGN model grid 
on the other hand actually tends to be too broad, 
demonstrating the sparsity of low metallicity AGN 
in the SP SS emission- l i ne ga laxy sample, as dis- 



cussed in Groves et al 



(|2006h . The low U NLR 
models are likely to be too faint to be observed 
in this sample. It should be remembered that the 
AGN and star forming galaxy models show are the 
extremes, of pure AGN contribution and pure H II 
regions. In between these models, seen in the right- 
hand wing of the diagram (the 'AGN' branch) are 
galaxies where both mechanisms, active star forma- 
tion and AGN, contribute to the galaxy spectrum. 

The second diagram (figure shows the mid- 
IR diagram of [Ne ni]15.6/im/[Ne n]12.8/im versus 
[S lv]10.5/zm/[Neli]12.8/mi. A surprisingly strong 
correlation between these ratios was found by 



Groves et al. (2008), who explored a large sample of 
emission line objects to see if the ground-observable 
[S lv]/[Nell] ratio could replace the space-only, ion- 
ization sensitive [Ne in] /[Neil] ratio. Here we show 



only a sub-sample of that collated bv lGroves et al 



(2008), displaying only star formation ionized re- 
gions (starbursts, blue compact dwarfs, extragalac- 
tic Hll regions, and galactic Hn regions), shown as 
crosses on figured Overlaying the observed data is 
same star-formation model grid of metallicity and 
ionization parameter as in figure[TJ as labelled in the 
colour bars. The model grid broadly covers the ob- 
servations, with both higher metallicity and lower 
ionization parameter towards the lower left of the 
grid (lower ionization). This is somewhat degener- 
ate with the age of the ionizing stellar population, 
with ages greater than ~ 4 Myrs moving strongly 
to the lower left of the grid. 

There are two interesting points to take from this 
grid. The first is that the models are unable to ex- 
plain the full spread of the data, though this may 
arise due to the ~ 2 0% uncertainty in the data (see 
IGroves et al. , 20081 for details) . The high values of 
[Siv]/[Neii] may come also from higher ionization 
parameters and or different abundance patterns, 
but the points to the right of the diagram suggest, 
as in the previous figure, that the low metallicity 
models (or youngest age or Wolf-Rayet models sim- 
ilarly) are lacking a hard enough spectrum to reach 
these ionization states. 

The second is that the observed relation between 
[Ne in] /[Neil] and [Siv]/[Neii] has a different slope 
and narrower spread than the model grid, suggest- 
ing that there may a higher order correlation going 
on between the metallicity and age of the stellar 
population and the ionization parameter of the ion- 
ized gas. 

The final diagram uses optical and far-IR lines 
of the same species, +2 , aiming towards the 
ionized gas diagnostics possible with Herschel 
Space Telescope, specifically the PACS spectrom- 
eter. Figure [3] shows the [Om] 52um/88/£nversus 
[O iiil 5007A/88/xm, based on the iDinerstein et al 
( 19851 ) diagram. This diagram is a strong diag- 



nostic diagram, separating out both density and 
temperature, while avoiding direct abundance and 
ionization effects. While unfortunately the den- 
sity sensitive 52^m/88/mi ratio cannot be observed 
in the local Universe with Herschel due to the 
lower wavelength limit, the temperature sensitive 
5007A/88/mi ratio can be obtained with the com- 
bination of Herschel PACS spectra with ground 
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Figure 3: Diagnostic diagram using only +2 emission lines, 
containing the density sensitive [Om] 52fj,m/88fim ratio ver- 
sus the temperature sensitive [O in]5007A/88/im ratio, that 
both use the Herschel observable [Oni]88^tm line. Shown 
are two different density (rtu = 10 cm -3 , lower left grid, and 
Tin = 350 c m" 3 , lower right gri d) instantaneous starburst 
models from Kewlev et alj (:2001), and two different density 
(nu = 0.1cm -3 , upper left green curve, and nn = 10cm -3 , 
upper right green curve) shock models from lAllen et alj 
(2008). In both models, the state of ionization (ionization 
parameter or shock velocity) increases from red to yellow. 
The four models are well separated in this diagram thanks 
to the diagnostic power of the ratios. 



based data. In reality however, caution must be 
taken with this ratio, as there are clear difficulties 
associated with the different wavelengths, such as 
matching the Herschel apertures and pixels with 
ground based data, and the large difference in dust 
attenuation between the two lines, both of which 
effects could overwhelm any small variations due 
to temperature. Nevertheless, to demonstrate the 
power of this diag ram we show t wo st ar forming 



model grids from iKewlev et alj (120011) . a nd two 



shock model curves from lAllen et all (2008) 



The iKewlev et ail (|200ll ) models are grids in 
metallicity (blue curves, same as Levesque mod- 
els) and ionization parameter (orange curves) us- 
ing a Starburst99 3 Myr old instantaneous burst 
as the ionizing source. What distinguishes the two 
model grids is the density of the ionized gas, with 
the left grid at nn = 10 cm -3 , and the right grid 
at nu = 350 cm -3 . The 52/^m/88/im ratio clearly 
separates the two grids. Though difficult to see due 
to their collapsed form, the 5007A/88/im ratio is 
affected by both the ionization parameter (weakly) 
and the metallicity (more strongly). However the 
diagnostic strength of this ratio is more clearly seen 
when compared to the much hotter ionized gas in 



(2008) model curves (green) 
show pure shock models at solar metallicity, mag- 
netic parameter of B/-y/n = 1/iG cm 3 / 2 , and veloc- 
ities ranging from 150 (red squares) - 1000 (yellow 
squares) kms -1 . Note that the velocities are shown 
as squares due to one of the features of ITERA, 
where if one of the parameters are set to a single 
value, in this case solar metallicity, the varying pa- 
rameter along the curve is marked by squares. As 
with the Kewley models, the two curves are again 
separated by density, with the left curve having a 
pre-shock density of nu = 0.1cm -3 , and the right 
nn = 10 cm -3 . What is clearly noticeable in fig- 
ure 02 is the clear vertical offset between the shock 
models and star forming models in the temperature 
sensitive [Om] 5007A/88/xm ratio, demonstrating 
clearly the possible capabilities of upcoming Her- 
schel-PACS spectral surveys. 

4. Summary 

We have presented here a new tool to enable as- 
tronomers to compare observations of emission line 
ratios with that determined by photoionization and 
shock models, ITERA, the IDL Tool for Emission- 
line Ratio Analysis. ITERA is an IDL widget 
tool which allows the user to plot ratios of any 
strong atomic and ionized emission lines as deter- 
mined by standard photoionization and shock mod- 
els. These line ratio diagrams can then be used to 
determine diagnostics for nebu lae excitation mecha- 



nisms (such as discussed in the lBaldwin et all 11981 



paper) or for nebulae parameters such as density, 
temperature, metallicity, etc. ITERA can also be 
used to determine line sensitivities to such parame- 
ters, compare observations with the models, or even 
estimate unobserved line fluxes. This tool, and as- 
sociated libraries and instructions, can be obtained 



at http : //www . brent groves . net/ it era . html 

At its core lies a library of emission-line nebulae 
models, c overing photoionization by young O and B 



stars (e.g. lKewlev et alll200lHDopita et all 12006) 



photoionizatio n by a hard power-law spectrum from 
an AGN (e.g. I Groves et all l2004bl) and ionization 



by shocks (e.g. lAllen et all 120081 ). These models 
broadly cover the space of observed emission-line 
objects. 

Finally, to demonstrate the capabilities of 
ITERA we present three different line di- 
agnostic diagrams; the classic BPT diagram 
of [Nn]AA6548,6584//ia versus [O iii]A5007/H/3, 
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showing the coverage of the star forming galaxy 
and AGN photoionization models when com- 
pared to SDSS emission-line galaxies, the mid- 
IR diagram of [Neiii]15.6/zm/[Nell]12.8/xm versus 
[Siv]10.5/Ltm/[Nen]12.8/im, revealing an interest- 
ing correlation of ionization parameter and metal- 
licity of stellar photoionization models when com- 
pared to observations, and the optical IR diagram of 
+2 , [Oni] 52/jm/88/im versus [Olll]5007A/88^m, 
demonstrating one of the diagnostic possibilities of 
Herschel. 

It is hoped that this tool will be of great use to the 
astronomical community when analyzing spectra of 
emission-line regions due to both its simplicity and 
ease of use, especially as the library expands with 
the growing number of photoionization and shock 
models becoming available. 
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